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The modified adiabatic scattering function previously used to compute 
phase shifts for positron- helium scattering below the positronium threshold 
is applied to the problem of annihilation in helium. A large enhancement factor 
relative to the Dirac rate is found, varying with energy from 2.30 to 3*1 6. 

The probability of finding the spectator electron in states of the He + ion 
other than the ground state is computed and found to be small but probably 


observable 
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I . Positron Scattering Function 

We have previously discussed the adiabatic approximation"^" and applied it 

p 

to the problem of low energy positron-helium scattering. Let us now consider 
its application to the computation of positron annihilation in helium. 

The scattering wave function for a positron of momentum k has the form: 


\(ri,r 2 ,x) = x k (x) [1 + G(ri,x) 








Here x is the positron coordinate, and £±>£ 2 . are the coordinates of the two 

electrons, all measured from the helium nucleus. The function 0 describes the 

ground state of the helium atom. The asymptotic form of the function Y will 

ik * x 

be correct if X v e “ “ since the correlation function G vanishes for large 

& 

values of x, and the two electrons are assumed to be independently polarized. 
Taking the positron interaction as perturbation 


V 



1 = 1,2 




( 2 ) 


one can evaluate G(r,^) correct to first order in V by the method of Dalgarno 

3 2 

and Lewis . If we make the shielding approximation employed previously 


0(r x ,r 2 ) = rr' 1 ^ 3 exp[-p (r x +r 2 )] 


( 3 ) 


the equation for G(r,x) becomes 


^ G . 2e g = £ ,[i . ^ ] - ase- 2 ^! + 


( 4 ) 


This is similar to the equation solved by Dalgarno and I^ynn for the case of 
hydrogen, and can be related to that result directly by making the change of 
variables y = |3r, z = |3x. Then one finds that 0G(y,z) satisfies the same 
equation as does the function derived by Dalgarno and Iynn , and can be taken 
over from their work. They derived an expression, in elliptical coordinates, 
which contains implicitly all terms in the Legendre polynomial expansion 

OO 

G(j£,z) = ^ gfiiy,*) P^cos a), [cos a = y*z]. ( 5 ) 

1=0 

Our previous experience^" indicated that the monopole term in the expansion 
gives excessive short-range attraction (correlation), and in the following it 
will be completely suppressed. Then, in elliptical coordinates one has the 
modified correlation function for helium 

0 (£,z) = p " 1 [f(z,X,p) - £ J d(cos 0 ) F(z,X,p)] = p _ 1 [F-F 0 ]. ( 6 ) 

-1 1 /o 

Here, X = z [y + (y 2 + z 2 - 2yz cos0) ' ], 1 ^ \ 

P = z " 1 [y - (y 2 + z 2 + 2yz cos0) ^ 2 ], -1 £ p s: 1, 

cos0 = (l + Xp) (X + p) 1 f 

F(z,X,p) = (X + p) A + [1 + 2Az -1 ] Xn(l + X) + D(z,p), ( 7 ) 

where D(z,p) = (l + z 1 ) [Ei(-z[l-p]) -in(l-p) 

-2z 

-e 


(Ei[z(l+p)] -4n[l+p])], 
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t" 1 dt. 


and A = ^ [ (l+z )e _2z -l] . 
II. Annihilation Rate 


The rate of annihilation of positrons in helium is given h y the expression^ 


1 



( 8 ) 


vhere r Q is the classical electron radius e 2 /mc 2 , N is the number density of 
atoms in the target gas, and is an effective electron number given by the 
expression 


\ = Jli d £i ^ <3^1 \(ri,r 2 ,x) f [6(ri-x) + 6(r 2 -x)3, 


(9) 


and provides a measure of the probahility that the positron and an electron 
are at the same point in space. If one were to set V = 0, thereby reducing 
the distortion of the atom as well as the positron scattering to zero, then 

... . ik«x , . 
Y k = 0 e *> (10) 


Thus the effective electron number = 2, and the annihilation rate goes over 
to the Dirac rate which is the quantum electrodynamic result for free particles. 
For our type of wave function. 


* \ 


dze Z I I 2 + G[z,z]) 2 + < G 2 >]^ 


vhere 


<Q 2 > = 


dri dr2 0 Q 0. 


(id 
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The energy independent parts of can be conveniently isolated by writing 
X k (x) in a partial wave expansion 

00 

X k (x) =£ x _1 U L (x) P l (cos 0)(2L + l)i L e i6L . (12) 

L=0 

Using 

00 

= Y j jj^kx) P L (cos0) (2L + l)i^ 

L=0 

one finds 

J <30 I X k (x) f = J dO | e 1 - ~ + (x k (x) - e 1 -**) f 

CO 

= 4tt[i + £ (2L + 1) (x' 2 U L 2 - j L 2 )], (13) 

L=0 

and 

00 03 

^ Zfc = 4 J dze~ 2z z 2 j^(l+G[ 2 ,z] -KG 2> J ^1+^ (2L+l)J^P 2 z 2 U 2 j 2 (kz/g 
o L=0 

(14) 


Thus, a rapidly converging series of energy- dependent corrections and a single 
energy independent term together yield the effective electron number. The 

2 

functions Ul are obtained numerically from the solutions discussed previously , 
and have the asymptotic form 


Uj^(x) -* k -1 sin (kx -In /2 + 6-^). 
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The problem thus requires the evaluation of <G 2 > and GCz,^) and the 
numerical integration of Eq. (l4). In elliptical coordinates, y = z is equiva- 
lent to X = n = 1, since y = z and cos0 = 1. For these values 


F(z,l,l) = 2A + [l + 2z -1 A]j&n2 

+(l+z 1 )[£hyz - e -2z (Ei[2z]-An2)], 


(15) 


where in y = 0.5771* Euler's constant. 

The integral in Eq. (6) takes the form 

i l 

^ Jd(cos6) F(z,X,|i)^ + ^ _ 2 = A Jd(cos0 ) 

-1 -l 

r 

+ [1 + 2z -1 A] J dX (X-l)An(l+X ) + J dki(l-|a)D(z,n) . (l6) 

i -l 

Then, after performing the integrations, one obtains the result 

G(z,z) = p' 1 {[l+z" 1 ][An 2yz - Ei(-2z) - l/2z + e' 2z (l+l/2z)] -lj- (17) 

which approaches zero as z -* 0, confirming the removal of the short-range 
parts of the correlation. Some numerical values of (l + G) are shown in 
Table I. 

The re ma i n ing term required to evaluate the annihilation rate is <G 2 > 
which has been evaluated numerically from the formula 
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<G 2 > = p" 2 [<F s > - <F^>] . (18) 

It also vanishes at z = 0 and decreases as a 4 for large z. For j3 = 1.5992 
it is also exhibited in Table I. In Fig. 1 and Table II the values of are shown over 

the energy range from zero to the positronium threshold. These have been com- 
puted numerically, using one set of solutions to the scattering [f3 =1.5992^^=0] 

p 

discussed previously , but including only the L = 0 part of the sum over partial 
waves. We estimate that the neglect of L > 0 contributes an error of less than 
2 $ to the annihilation rate, and the error vanishes for k -* 0. The rates 
obtained are very high , and exhibit considerable variation with energy, which 

7 

should soon be observable experimentally. 

III. Excitation of the Residual Ion 

Suppose that an electron is annihilated when the positron is at position 

x.' The annihilation takes place very rapidly compared with atomic times, and 

so the subsequent behavior of the residual He ion can be treated by the sudden 

8 + 

approximation . The probability that the He ion will be in the state whose 
wave function is 0 n is then given by 

P n = l^l a n ( ^ f, ^9) 

n 

where a^(x) = Jdr 0*(r) \(x,x,r) . 


The denominator in Eq. (19) is 
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JJJ d £ d £’ ^ K { l"> 0 S ( ^ } ] 

•* n ** 


( 20 ) 


and since {0 n } is a complete set of one-electron functions, closure can be 
applied to the sum in Eq. (20). Thus 


and 


P 

n 


J I *n(x) f /J J te&r | \(x,x,r) f 
Jdx I an(x) ?/ 



= 1 . 


( 21 ) 


If the probability for exciting states other than the ground state is 
not too small, there may be a possibility of observing the radiative de- 
excitation of the residual He + ion, following annihilation. A favorable 

O 

transition might be the 2p - Is ultraviolet line at 304 A, and we now proceed 
to compute its intensity, using the modified adiabatic wave function. 

From Eqs. (l) and (19) one obtains 

p n = 2(TTZ k ) _1 Jdz e -2z | x k f [(l + G[z,z])Q n + w n ( z )] 

where On = (P 3 /tt) ^ J dr e" Pr 0 n (r) (22) 

W n (z) = 0 3 /tt ) / 2 J dr G(y,z) 0 n (r). 

For excitation to the 2p state of the ion = 0 and only the dipole part 
of G contributes to W n (z). This is obtained by inserting Eq. ( 5 ) into Eq. (4) 


and yields 
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gi(y,z) = P" 1 [z 2 [y + y 2 ^] + C x [e'^y 2 + (l + y _1 + ^y" 2 )] } 


for y < z 


= P _1 {z+ C 2 [l + y" 1 + ^y' 2 ] j- 


(23) 


• - 


for y > z 


where Ci = § e~ 2z (l + z” 1 ) 2 , 

C 2 = Ci - z + § (l - z" 2 ). 


and then 

00 

W 2 p(z)=33 ^ 2 J ayy 3 e"^ 1+P ^ y gi(y,z). (24) 

o 

The probability Pg has been computed for the same parameters as before and 
is found to be essentially independent of positron energy, (for L = 0 scatter- 
ing) and equal to 0.82$ . Thus at least 0.82$ of all annihilations should be 
followed by the 2p -* Is transition. In addition, however, the experimentally 
observable radiation should include the effect of cascades from other states 
which pass through the 2p state. S-states should dominate the cascade process. 

For S-states, W n = 0 since no monopole term is contained in G. The S- 
state probabilities are 

P ns * t 1 + R ] 1 


where R = J d:z e' 2z | f <G 2 >/Jd:z e~ 2z | X k f [l + G ] 2 


( 25 ) 


Invoking closure over S- states, one finds that P(all S) = [l+R] -1 = 99 $ where 
R has been numerically evaluated, and proves to he essentially energy- indepen- 
dent . Thus 99 $ of all annihilations give S- states of the ion, and 0.82“$ go 
to the 2p state, leaving a negligible fraction in the higher p or L > 1 states. 

The partition among the various discrete S-states is then given by Q n a . A 

9 

general formula due to Mittleman is 

0/ = (8p/n) 3 (2-p) 2 (P-2/n) 2n-4 (p+2/n)" (32n+4) , (2 6) 

which was used to prepare Table m, which also shows the branching ratios^ 
to the 2p state. One finally finds that 1.2$ of all annihilations should be 
followed by the 2p Is transition. In addition, if a sufficiently strong 
electric field is applied to the target gas, the metastable 2s level will mix 
with the 2p level 11 , and an increase in the 2p Is transition rate to 3«2$ of 
the annihilation rate will occur, and should be observable. 

IV* Discussion 

We have presented above a calculation of the energy-dependent annihilation 
rate of positrons incident on a helium target, as well as a prediction of the 
intensity of the He + ultraviolet radiation which is expected to follow the 
annihilation. In this work a simple model wave function describes both pro- 
cesses, as well as the elastic scattering. The predictions are thus not 
expected to be exact, but previous experience has given some cause for optimism. 

Nevertheless the results derived here must be used cautiously, as a guide 
only, until such time as more accurate solutions of the positron-helium scatter- 
ing problem become available. As an example, we may note that the value of 



the parameter 3 = I.5992 used throughout has "been chosen to produce agreement 

-4 12 

with the long-range x potential , and thus does not necessarily represent 
well some of the other properties of the system. In a somewhat different 

Q 

context, Mittleman has used a five-term exponential approximation to describe 
the helium atom wave function, and finds that a more accurate value of the 2s 
excitation probability is 2 .23%, rather than our 2 . 00 $ . This agreement indi- 
cates that our qualitative results should be used as preliminary estimates in 
designing experiments, which may eventually test the future detailed theories. 
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Table I. Numerical values of correlation functions required for the calculation 
of the positron- helium annihilation rate, for full monopole suppression 
and p = 1.5992. 


z [1 + G(z,z)] 2 <q 2 > 


0 

u — « \ ~ J <-• / J 

1 

0 

0.5 

1.886 

.0230 

1.0 

2.733 

• 0473 

1-5 

3.447 

.0495 

2.0 

4.028 

.0397 

2-5 

4.503 

.0282 

3.0 

4.900 

.0188 

3.5 

5-241 

.0123 

4.0 

5.541 

.0080 

4.5 

5.810 

.0053 

5-0 

6.055 

.0036 


0 


Table II. Effective electron number for annihilation at selected values of 
incident positron energy. The Dirac rate corresponds to = 1. 


Energy (eV) 


0 



3-16 

0.54 

2.79 

1.22 

2.58 

2.17 

2.43 

3.40 

2.35 

4.90 

2.31 

6 .66 

2.30 

8.70 

2.33 

11.0 

2.38 

13-6 

2.44 

17.8 

2.54 


Table IEE- Population and branching ratios following annihilation. The pro- 
bability for exciting the various S-states of He + listed in column 
1 are shown in column 2, and the fraction: of each which eventually 
reaches the 2p state is shown in column 3* [Cascades like 5 s 4p “* 

3s -* 2p are included, but are very small] . Column 4 gives the 
observable population of the 2p state due to each S- state. 

Contribution to 


State Population Branching Ratio to 2p 2p Population 


Is 

95*34$ 

1 

0 

0 

2s 

2.00$ 

0 

0 

3s 

0.30$ 

1 

• 30$ 

4s 

0.10$ 

0.58 

.058$ 

5s 

0.04$$ 

0.45 

.022$ 

6s 

0.027$ 

0.43 

.012$ 




Total = .39$ 


Figure Captions 


Effective electron number for 
energy in eV. The Dirac rate 


annihilation versus incident positron 
corresponds to 1 / 2 Z, = 1. 
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